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Abstract 

We discuss the ICARUS T600 detector capabilities in electromagnetic shower 
reconstruction through the analysis of a sample of 212 events, coming from the 
2001 Pavia surface test run, of hadronic interactions leading to the production of 
TT^ mesons. Methods of shower energy and shower direction measurements were 
developed and the invariant mass of the photon pairs was reconstructed. The (7,7) 
invariant mass was found to be consistent with the value of the 7r° mass. The 
resolution of the reconstructed vr*^ mass was found to be equal to 27.3%. An improved 
analysis, carried out in order to clean the full event sample from the events measured 
in the crowded environment, mostly due to the trigger conditions, gave a mass 
resolution of 16.1%, significantly better than the one evaluated for the full event 
sample. The trigger requirement of the coincidence of at least four photomultiplier 
signals favored the selection of events with a strong pile up of cosmic ray tracks and 
interactions. Hence a number of candidate tt^ events were heavily contaminated 
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by other tracks and had to be rejected. Monte Carlo simulations of events with 
vr'^ production in hadronic and neutrino interactions confirmed the validity of the 
shower energy and shower direction reconstruction methods applied to the real data. 

Key words: Decay of tt mesons, Calorimeters 
PACS: 13.20.Cz, 29.40.Vj 



1 Introduction 

The Liquid Argon (LAr) Time Projection Chamber (TPC) is an option for 
the next generation of very large mass neutrino detectors. The LAr TPC was 
proposed by C.Rubbia [1] in 1977 and, after many years of extensive R&D 
studies the T600 ICARUS detector [2] has been constructed. The read-out of 
ionization electrons traveling in a uniform electric field provides the spacial 
reconstruction of an event and allows for calorimetric measurement of the 
energy release. In year 2001, the T600 ICARUS detector was extensively tested 
in Pavia (Italy), with an exposure to cosmic rays at the Earth surface. During 
over 100 days of continuous data taking, about 30000 cosmic ray events have 
been collected. In a series of papers the ICARUS Collaboration presented the 
following results from the T600 detector test run: observation of long ionizing 
tracks [3], LAr purity [4], electron recombination in LAr [5] and measurement 
of through-going particle momentum by means of multiple scattering [6] . The 
analysis of the sample of stopping n events led to the measurement of the 
decay spectrum, and allowed to determine the energy resolution of the T600 
ICARUS detector for the measurement of electrons with energy lower than 
about 50 MeV [7]. 

The ICARUS T600 detector is at present in the final phases of re-assembling in 
the underground Hall B at INFN Laboratori Nazionali del Gran Sasso (LNGS) , 
Italy, and its commissioning is foreseen to start within the year 2008. 
In this paper we discuss the ICARUS detector capabilities in electromagnetic 
shower reconstruction through the analysis of a sample of events with hadronic 
interactions leading to the production of 7r° mesons. The decay of a particle 
with known mass gives a good test for the energy calibration. Another 
motivation for such studies is the following: in order to distinguish charge 
current (CC) interaction 

Ve + n ^ e~ +p, (1) 
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from the neutral current interaction 

Vn + n ^ + + hadrons, (2) 

a very good detector capability to separate electrons and pions is required. 
The tt" decays in about 98.8% into two photons. The Dalitz decays, leading to 
the production of one photon and electron-positron pair, compose about 1.2%. 
Both photons and electrons may create electromagnetic showers and therefore 
lead to the similar event topology for both NC and CC interactions. 
The paper is organized as follows. The ICARUS T600 detector set-up is briefly 
described in Section 2. Section 3 concerns the search and the selection of 
cosmic ray hadronic interactions in T600 detector with the production of the 
7r° meson. In Section 4 the methods used for the measurement of the energy 
released in LAr by the two electromagnetic showers coming from 7r° — > 7 7 
decay and for the reconstruction of tt" mass are illustrated. The analysis of 
the selected event sample from real data is presented in Section 5. In Section 

6 the precision in the measurement of the electromagnetic shower energy 
deposition in LAr is evaluated in the condition of the test experiment through a 
simulation of electrons and photons showering inside T600 detector. In Section 

7 the Monte Carlo simulations of vr'' in ICARUS T600, used to validate the 
reconstruction routines for electromagnetic showers, are shown. In Section 8 
the same method of event reconstruction used for the real data is applied to 
Monte Carlo simulated muon neutrino interactions with 7r° production. The 
possible mixing of such events with the electron neutrino CC interaction is 
discussed. The work is summarized and concluded in Section 9. 



2 Experimental set-up 

The T600 ICARUS detector is composed of two identical adjacent "T300" half- 
modules filled with Liquid Argon. Two Time Projection Chambers (TPCs) 
with common cathode placed in the middle, the electric field shaping system, 
LAr purity monitors and photo-multiplier tubes (PMTs) are enclosed in each 
half-module. The external thermal insulation, together with a Liquid Nitrogen 
circuit, allows to stabilize the LAr temperature. The required LAr purity is 
achieved by a purification system. Each T300 half-module has the following 
internal dimensions: 3.6 m (width), 3.9 m (height) and 19.9 m (length). Along 
the longest side walls (left and right) of each half-module, three parallel planes 
(two Induction planes and one Collection plane) of anode wires are placed. 
They are parallel to the cathode and form two TPCs. They are oriented at 
60° with respect to each other, 3 mm apart, with wire pitch of 3 mm. The total 
number of wires in the T600 detector is 53248, whereas the volume of LAr in 
the instrumented part of the detector is 340.3 m^ (476.5 tons). The ionization 
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electrons, forced by a uniform electric field of 500 V/cm perpendicular to the 
cathode, drift toward the anode wires inducing a signal on those which they 
pass by. The appropriate voltage biasing causes non-destructive signal to be 
induced first on two Induction wire planes, whereas the ionization electrons 
charge is finally collected in the Collection wire plane. By combining the 
information from two-dimensional projections, provided by each pair of wire 
planes, with the position along the drift time, it is possible to obtain a three- 
dimensional reconstruction of an event. The zero time to of an event is defined 
by the PMT detection of LAr scintillation light. The position along the drift 
time is given with respect to tg. 

Several triggers have been applied during the T600 ICARUS detector test run. 
An external trigger system made of plastic scintillators was used for detection 
of muon tracks. An independent trigger system was based on photo- multipliers, 
which detect scintillation light produced during interaction of charged particles 
with LAr. Most of the candidates used in this analysis were found in 
the data collected with this internal trigger. The detailed description of the 
T600 ICARUS detector design, construction, and performance is given in [2], 
where the conditions of the technical run and the event off-line reconstruction 
procedure are discussed as well. 



3 Data selection 

Interactions in LAr producing 7r°s decaying into photons have been selected 
through a visual scanning for events with at least two showers pointing to 
the interaction vertex: Fig. 1 shows an example of an ICARUS T600 full view 
event, with a tt*^ candidate that was selected during the scanning phase. Most of 
the selected interactions are initiated by hadrons, hence we also use the term 
"hadronic interactions". Altogether, about 7500 images have been scanned. 
The majority of the events were collected with triggers based on signals from 
photo-multipliers, rich in hadronic interactions. 

It has to be stressed that the T600 ICARUS test run internal trigger based 
on photomultipliers was aimed at air showers, so that many events were not 
measurable because of their complexity due to a crowded environment. The 
surface exposure has in fact as a result that a number of candidate 7r° events 
can lie in zones of the detector of strong pile up of cosmic ray events, which can 
partially overlap the electromagnetic showers of the photons of the tt" decay. 
This affects either the capabihty of the visual scanning in disentangling the 
good events from an environment populated by many through-going particles 
or the precision in the shower energy reconstruction for selected events with 
partial overlap with cosmic ray events. 

The following scanning rules were adopted for the pre-selection of hadronic 
interactions with tt^ production: 
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• At least two, well separated electromagnetic showers originated by photons 
pointing to the interaction vertex are demanded. Shower verteces and 
directions are defined by the initial electron-positron pairs, so that it is often 
easy to decide if they are pointing to the interaction vertex. The showers, as 
well as the interaction vertex, must be clearly visible in at least two views. 

• Knowledge of the event zero time (to) is mandatory in order to determine 
the interaction vertex and the shower spatial positions in the chamber, as 
well as to apply proper corrections to the collected signals for the charge 
attenuation due to the finite electron lifetime in Argon, associated to the 
electron attachment process to electro-negative impurities. 

A unique to determination is possible in the following cases: (1) at least one 
track (but not showers) from the interaction is crossing the anode wires {to — 
fjLS for crossing points); (2) at least one track (but not showers) is crossing 
the cathode plane {to — 962 /is for crossing points); (3) an unique association 
between the interaction and the PMT induces signals on Collection wires 



• Interactions with neutral meson candidates are rejected when ionization 
due to other tracks and interactions in their neighborhood is so high that 
it seriously infiuences the shower energy determination. 

Every candidate event was analyzed independently by (at least) two groups 
from different laboratories in the Collaboration. After all these requirements 
altogether 212 hadronic interactions with at least one candidate neutral meson 
were collected for further analysis. As said before, the crowded environment 
of ICARUS T600 test run data reduced the usable data sample substantially: 
examples of accepted and rejected events with neutral meson candidates are 
shown in Fig. 3 and in Fig. 4, respectively. 



4 Data reconstruction 

The invariant mass M of a particle decaying into two photons is given by 



where Ei and E2 are the energies of the two photons and 612 is the angle 
in between. In our case, the goal was to measure the energies of two 
electromagnetic showers coming from the photons and to reconstruct the 
photon directions. 

In Fig. 5, one of the 212 selected events is shown: in both Collection and 
Induction views an hadronic interaction vertex and the two shower starting 



(Fig. 2). 




(3) 
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Fig. 1. Example of an ICARUS T600 full view event (Run 939 Event 130). On the 
left side, close to the center of the image, a large hadronic interaction leading to the 
production of a number of secondaries is visible. Along all the image m.i.p. tracks 
can be seen. A zoomed view of an hadronic interaction leading to vr'' production and 
subsequent decay in two photons is shown. The scale of the image does not match 
the real one. 

points are visible. The drift position coordinate T, where T(cm) = drift time 
(s) X electron drift velocity (cm/ s), is common for the two 2D projections of the 
events, while the U (cm) and V{cm) wire coordinates indicate the wire position 
in Collection and Induction views, respectively. The U, V and T coordinates 
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Left collection view 
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Fig. 2. An example of event in which the to is uniquely determined from PMT 
induced signal on Collection wires, visible at the drift coordinate T~ cm as an 
ionization strip along the horizontal U wire coordinate (Run 700 Event 58). 



Left collection view , , Right collection view 
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Fig. 3. Examples of events from the final data sample: a) ir^ candidate (Run 939 
Event 130); b) a multi-shower event (Run 939 Event 190). 

of the interaction vertex and those of the e.m. shower starting points are used 
for the three dimensional reconstruction of the photon direction, hence of the 
angle 9i2 between the two photons of the 7r° decay. 

The basic step for the calorimetric reconstruction of electromagnetic showers 
is the hit finding and reconstruction procedure. We define hit as the portion 
of particle track which energy is read by a given wire. The hit search and 
reconstruction is performed on the raw signal of the Collection wire plane 
following the standard procedure for ICARUS data analysis described in [2]. 
After hit reconstruction, the following parameters are available for each hit: 
its wire-drift position U and T and its area (ADC x jis). 
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Fig. 4. Examples of candidates not included in the final data sample: a) neutral 
meson candidate with overlapping photon showers(Run 975 Event 65); b) candidate 
event in a crowded environment (Run 710 Event 29). 

A dedicated HIGZ-based Graphical User Interface visualization tool allows to 
select directly the region with the hits belonging to a shower, by building a 
polygonal line around the shower directly on the Collection view by means of 
interactive calls of HIGZ ambient. Then hit finding routine is applied within 
the bounds defined by the polygonal line (sec Fig. 6). It has to be noticed that 
the definition of the area surrounding the shower is left to the user choice, 
so that small variations of the reconstructed energy can be found between 
measurements done by different users. This happens especially in the case the 
two photons of the 7r° decay are emitted at small angle, so it can be difficult 
to decide whether an energy deposition belong to the first shower or to the 
second one. 

Once all hits belonging to a shower are found, the overall shower energy E is 
computed by summing all energies where i is the index of the hits inside 
the shower area: 



1 CW 

E,{MeV)^—--A-Q,, (4) 

£ K 

where: C = (152 ± 2) x IQ-^ fC/(ADC x /^s) is the calibration factor [4]; 
W — 23.6^0 1 eV is the average energy needed for the creation of an electron- 
ion pair [8]; = 0.640 ±0.013 is the electron- ion recombination factor [7]; A = 

e -^e is the drift electron lifetime factor: tj is the hit drift time coordinate and 
to is the event zero-time, is the actual electron lifetime in LAr, that varies 
run by run [4]. Its value as a function of the acquisition date was monitored 
during the T600 data taking [9], with relative errors at the level of 5%; finally, 
Qi is the hit area. The shower energy E is then corrected by multiplying it 
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Fig. 5. Example of selected candidate vr*^ — > 7 7 event (Run 712 event 7). Hadronic 
interaction vertex and e.m. shower starting points are indicated in both Collection 
and Induction views. 

by a factor coming from Monte Carlo simulations of e.m. showers from 
photons and electrons in the range 50 ^ 5000 MeV in ICARUS T600, taking 
into account the shower energy reconstruction efficiency (see Sec. 6). It has 
to be noticed that the correction for the reconstruction efficiency in eq. (4) 

depends on the attenuation factor A = e ~ coming from the drift electron 
finite lifetime of the Monte Carlo generated events, so that the e value changes 
shower by shower, ranging from e = 0.69 ± 0.03 {A = 0.5) to e = 0.83 ± 0.04 
{A = 1) (see Fig. 11). The error on the shower energy is then evaluated from 
the propagation of all contributions in eq. (4). 

The angle 612 between photons is then evaluated from the spatial coordinates 
of hadronic interaction vertex and of starting points of the showers as explained 
above. 



5 Results 



5.1 (7, 'j) invariant mass 



Most of the 212 events forming the sample were measured in several 
laboratories. The difficult condition of data analysis due to the presence of a 
crowded environment sometimes has as a result that the same event measured 
in different laboratories produces different values in vertex coordinate and 
shower energy determination. As explained above, the somehow arbitrary 
choice of the user to build the polygonal line surrounding the e.m. photon 
shower of the 7r° decay could in principle introduce some systematics in the 
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Left collection view Left Induction II view 
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Fig. 6. Example of the application of the routines that allow to select the hits 
belonging to an e.m. shower: the polygonal lines are built by the user. Reconstructed 
hits inside the areas defined by the polygonal lines are marked with red dots (Run 
712 event 7). 

(7,7) invariant mass distributions coming from different laboratories. Fig. 7 
shows the (7,7) invariant mass distribution of three sub-samples of events, 
each provided by a different laboratory involved in the data analysis. It appears 
that on average no systematic differences between measurements carried out in 
different laboratories are present. At this stage of the analysis, the e parameter 
taking into account the shower energy reconstruction efficiency, see eq. (4), 
has not been considered yet, hence the mass values are underestimated with 
respect to the nominal ones. As a remark, it has to be noticed that in LNGS 
underground condition the events will be very clean, so that each 7r° decay 
event will be measured just once unambiguously. The event sample for the 
following analysis was then built by taking just one measurement per event 
chosen randomly from the available ones in the case of events measured by 
more than one laboratory, then adding to the sample all the other events 
measured just once. 

A preliminary consistency check was performed on this event sample: the 
distribution of distances of the shower starting points from the corresponding 
interaction vertex is related to the radiation length Xq in LAr, Xq = 0.14 m 
[11]. This is shown in Fig. 8. The value of the mean distance obtained from a 
fit with an exponential function is: 



= (0.174 ±0.008) m, (5) 

that is in agreement with the expected length for photons X^ = | ■ Xq = 0.18 
m for LAr. 

A fiducial volume cut was then applied by requiring that at least the first 
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three radiation lengths of both showers were fully contained in the chamber. 
This reduced the event number to 164, of which 126 were two-shower and 38 
were multi-shower events; the final number of vr*^ candidates in the selected 
event sample was then found to be equal to 196. The invariant (7,7) mass 
distribution of this 196 candidate sample is shown in Fig. 9, as well as the 
7r° energy distribution. The correct mass combinations of multi-shower events 
were disentangled after the minimization of the parameter: 



X 



i=l 



m. 



i\2 



(6) 



evaluated for all the possible combinations of showers. (771,^0 = 135 
MeV/c^ [11])- Here, n is the maximum number of vr^s allowed (1 for a three 
shower event, 2 for a four and five shower event, 3 for a 6 shower event). The 
average mass given by a gaussian fit is: 



■77 



134.4 ± S.OMeV/c^ 



(7) 



The evaluation of all systematic contributions to the error on the vr" mass (see 
Sec. 4) introduces an additional contribution to the mass error of 7.1% (9.5 
MeV/c^). The RMS uncertainty of the reconstructed mass was found to be 
equal to a = 36.7 ± 3.4 MeV/c^ (27.3%). Thus the result, in agreement with 
the tt" mass, demonstrates that our mass reconstruction method is reliable. 



400 500 600 

mass(MeV/c') 



Fig. 7. (7,7) invariant mass distributions (solid, dotted and dashed lines) of three 
different laboratories involved in the data analysis. 



13 




Range (cm) 



Fig. 8. Distribution of the distances between shower starting point and interaction 
vertex for LAr radiation length evaluation. 
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Fig. 9. (7,7) invariant mass (a) and energy (b) distributions for the sample of 196 
7r° candidates. The average mass is m^^ = 134 A ± 3.0 MeV/c^ {a = 36.7 ± 3.4 
MeV/c^). The mean value of the 7r° energy is E = 707.7 MeV. 

5.2 Crowded environment suppression 



Although on average the measurements done in different laboratories are 
in agreement, see Fig. 7, nonetheless for this particular analysis of events 
collected on surface by exposure by cosmic rays it has been chosen to improve 
the results by cleaning the 212 event sample from the events measured in 
a crowded environment. To do so, further rejection cuts were related to 
the consistency of the two independent measurements performed in different 
laboratories. The purpose of this refined analysis is to show that in an 
environment more similar to the one expected in the underground laboratory 
the 7r° mass distribution has to be narrower than the one found for the full 
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sample. 

The interaction vertex and shower starting point coordinates are then 
considered to be in agreement if the absolute value of the differences between 
two independent measurements of all three coordinates is smaller than five 
times the width of the distribution of such differences. Then, if the relative 
difference %\+e2) between measurements of the shower energies is smaller than 

2 

50%, then the measurements are considered to be in agreement, so that only 
these are used to built the final mass distribution. The values used for these 
cuts largely suppress the 7r° decay events belonging to a crowded environment 
while keeping a relatively high statistics for the following analysis. For events 
which passed these tests the average values of the measured positions and 
energies were calculated and used for 7r° mass determination. 
The sample which passed these selection criteria was then composed of 89 
events with at least one 7r° candidate; 72 were two-shower and 17 were 
multi-shower events. As done for the global data sample, the correct mass 
combinations of multi-shower events were disentangled after the minimization 
of the parameter, see eq. (6), evaluated for all the possible combinations of 
showers. After this procedure, a final sample of 97 7r° candidates was available 
for the analysis. 

The (7,7) invariant mass distribution for the final sample of 97 7r° candidates 
is shown in Fig. 10a together with the total 7r° energy in Fig. 10b. The average 
mass given by a gaussian fit was: 

m^^ = 139.9 ± 2.8MeV/c^ , (8) 

with a systematic contribution of 9.9 MeV/c^. The 7r° mass resolution is now 
equal to 16.1% (a = 22.6 ± 2.9 MeV/c^), that is 8.9% less than the one 
evaluated for the sample without consistency cuts. This demonstrates that in 
the particular case of the ICARUS T600 Earth surface data, a refined analysis 
was necessary in order to take into account the strong contribution of the 
crowded environment surrounding some of the 7r° production and decay events. 
These crowded events are badly measured and bring to major differences 
between measurements carried out in different laboratories. 



6 Energy resolution for electromagnetic showers 

In this section an analysis of the energy resolution for electromagnetic showers 

in the energy range 50 -i- 5000 MeV is presented. The energy of Monte 
Carlo simulated electrons and photons inside the ICARUS T600 detector was 
reconstructed with the same tools used for the 7r° events. The simulation 
was performed using the FLUKA [10] program for the samples of electrons 



15 




M^ii {McV/c' ) Ene rgy (MeV) 



Fig. 10. (7,7) invariant mass (a) and energy (b) distributions for the final sample 
of 97 vr'^ candidates obtained after crowded environment suppression. The average 
mass is m^^ = 139.9 ± 2.8 MeV/c^ [a = 22.6 ± 2.9 MeV/c^). The mean value of 
the 7r° energy is £■ = 641.8 MeV. 

and photons generated at six different energies with different orientation 
in the detector. The measured average energies appeared to be shifted 
with respect to the true energies defined at the generation. To obtain an 
agreement between generated and reconstructed energies, the shower energy 
reconstruction efficiency (i.e. the ratio between measured and generated 
shower energy) e has to be taken into account. The inefficiency can be due 
to misidentification of the hits by the hit finding algorithms, as well as to 
the loss of the tail of the reconstructed tracks, whose energy deposition can 
become lower than the intrinsic wire background noise, especially in the case 
of low values of the electron lifetime Xg: in this case the drift electron finite 
lifetime factor A = e ~ is crucial. In Fig. 11 the reconstruction efficiency 
e as a function of A is shown, measured for a sample of Monte Carlo events 
at several values of Te'. it is found that e rises with A, and that it is well fitted 
with a fourth degree polynomial. Despite the arbitrariness of the choice of 
this particular function for the fit of e versus A, it was found that choosing 
other types of curves for the fit did not introduce systematics in the value of 
e, the fluctuations of which remain always well below the 1% in the region A 
= (0.5-i-l.O) where the real data lie. 

Depending on the values of this attenuation factor measured in the real data, 
e was found to range from 0.69 ± 0.03 {A = 0.5) to 0.83 ± 0.04 {A = 1). 
Moreover, in the range of e.m. shower energy analyzed, no dependence of e 
from the shower energy, the attenuation factor being the same, was found. The 
effect of the e parameter taking into account the shower energy reconstruction 
efficiency is shown in Fig. 12, where the dependence of the reconstructed vr'' 
mass from the distance between the main hadronic interaction vertex and the 
wires is shown (sample of the 196 7r° candidates). When the e efficiency is 
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not taken into account (red dots), an effect of underestimation of the mass, 
especially near the cathode {y = 150 cm), is evident. After the shower energy 
correction (green dots), this effect is much smaller and, as expected, the 
correction is larger for events far from wires {y > 60 cm). This effect can 
be seen in a more quantitative way by giving the values of the x^s from 
the fit of the two distributions with a flat line. The value of x^/^df after 
the shower energy correction {x^/ndf = 1.55/4) is significantly smaller than 
the one without correction {x^/ndf = 5.52/4), demonstrating that after the 
correction the mass distribution is almost flat. 




Attenuation factor A 

Fig. 11. Reconstruction efficiency e as a function of drift electron finite lifetime 
factor A: data have been fitted with a fourth degree polynomial. 
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Fig. 12. Dependence of the reconstructed vr'' mass from the distance between the 
main hadronic interaction vertex and the wires with (green dots) and without (red 
dots) the correction for the e reconstruction efficiency parameter. The red (green) 
straight represents the result of a fit of the red (green) dots with a flat line. 

The energy resolution at fixed energies is given by the RMS of the 
measurements. The following formula was fitted to the experimental points: 
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(9) 



where PI and P2 are the free parameters of the fit. The first part of the formula 
represents the typical electromagnetic calorimeter behavior. The introduction 
of P2 parameter improves the fit quality. Its presence can be due to the fact 
that the amount of low energy signals at the end of the showering process grows 
with the initial energy. These signals can be mixed with the overall background 
when the hits belonging to the shower are selected. Fig. 13 illustrates the 
resulting energy resolution of the T600 ICARUS detector in the reconstruction 
of the electromagnetic showers. 
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Fig. 13. as a function of V^. Dashed line - fit to eq. (9) (P2 = assumed), 
solid line - fit to eq. (9). 

The best fit with the formula of eq. (9) was obtained for the following 
parameter values: PI = (0.33 ± 0.03) VMeV and P2 = 0.012 ± 0.002. 
The energy resolution in LAr detector was previously studied on a base of 
a sample of MC generated low energy electrons (E < 50MeV) [7]: the best fit 
was obtained for PI = (0.11 ± 0.01) VMeV and P2 = 0.025 ± 0.03. However 
in that case PI and P2 were correlated because of the small energy interval 
involved. The worse resolution for the high energy electromagnetic showers 
can be explained by larger fluctuations of integrated background because of 
the larger area covered by the shower as well as to the hits overlapping in 
the regions of large density of electron tracks. Moreover, in the case of a high 
energy shower it is not obvious which surrounding signals still belong to the 
shower. The energy resolution has been obtained for the specific method of 
energy reconstruction based on calorimetric information from Collection view 
only. 
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As a check of the rehabihty of the energy reconstruction method, the 
distribution of the so called prills : pull{m^a) = ^"^^^^ ""j" ^ was calculated 
for the real data, using eq. (9) to calculate the error (T{m^^) as a function 
of a{E^). The RMS of the distribution was found to be equal to one within 
errors, confirming the validity of the formula obtained on the base of the 
measurements of the photon showers. 



7 Monte Carlo simulation of tt production from hadronic 
interactions 

Monte Carlo simulations of three different samples of 7r°s with different 
kinematic and event environment conditions were carried out. The events 
have been generated using the FLUKA full simulation of T600 detector. 
They allowed to confirm the vahdity of the 7r° mass reconstruction tools in 
environments similar to the one coming from real data. The three generated 
7r° samples were: 

(1) 7r° decays at rest; the invariant mass is reconstructed as a sum of the 
photon energies since the relative angle is exactly equal to 180°. Only 
errors on the photon energies contribute to the mass resolution. 

(2) E"*" — > 7r° + p with subsequent tt'' decay; the pion from the two body 
decay at rest has fixed momentum of 189 MeV/c and the vertex of 7r° 
decay is precisely determined by the position of the low energy proton 
from S"*" decay. Shower verteces are also well determined because no 
tracks overlap and affect their recognition. 

(3) 7r° decays observed in interactions of 7r~ mesons with Argon nuclei 
at 2 GeV incident energy. The events resemble most our real events, 
corresponding to hadronic interaction. 

Examples of each class of these events are shown in Fig. 14. The (7,7) invariant 
mass distribution for the three measured Monte Carlo samples is shown in 
Fig. 15a, 15b and 15c. It can be seen that the position of the m-y-y peak is 
consistent with the pion mass for all three generated vr" samples. 
The RMS of the invariant mass distributions presented in Fig. 15 are equal to 
(7.4 ± 1.3) MeV, (10.7 ± 4.8) MeV and (25.8 ± 4.7) MeV. The evident increase 
in RMS is mainly due to the event crowded environment, which augments the 
difficulty of the event interpretation because of the higher energy and the 
consequent presence of other reaction products. Moreover, in the decay at 
rest the relative angle between photons was not measured, but assumed to be 
equal to 180°. It has to be noticed that Monte Carlo events were generated 
in the center of the detector. For real events the effect of some energy leakage 
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Fig. 14. Examples of neutral pion decays as seen by FLUKA-ICARUS detector: 
decay at rest (a), tt^ from decay (b) and from vr" interaction in LAr (c). 

through the walls of the detector can appear despite of the fiducial volume 
cut. 
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Fig. 15. (a) m^^ for pion decays at rest: m^^ = 133.3 MeV/c^ (RMS = 7.4 MeV/c^ 
(b) m^^ from S+ — > 7r^+p decay: m^^ = 137.8 MeV/c^ (RMS = 10.7 MeV/c^j; {. 



m^j from 2 GeV vr LAr interactions: m^y 



130.6 MeV/c^ (RMS = 25.8 MeV/c^ 



8 Simulation of tt^ production in u interactions 



In the previous sections it has been shown that the results of the analysis of 
real and simulated events with 7r° production are consistent. Therefore, the 
same method of generation, measurement and reconstruction was applied to 
simulated neutrino-induced pion production to check the potentiality of the 
LAr technique in this respect. An important point in the future experiments 
looking for ^13 is the efficiency of distinction between photon-induced showers 
from 7r° decay and the electron showers from electron neutrino interactions. 
The neutral pion production from NC neutrino interactions in the Liquid 
Argon ICARUS detector was simulated using NUX-FLUKA chain of 
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programs. The neutrino energies were chosen in a region of maximal oscillation 
probability between 1 and 5 GeV. 

The measured distribution of the two photons invariant mass is presented in 
Fig. 16. To disentangle the proper photon pair combination for events with 
four and six photons pointing to the same vertex, the pairs with minimum 
of the difference between measured values of (7,7) invariant mass and the 7r° 
mass were selected, see eq. (6). The position of the maximum was found to 
be consistent with the vr^ mass value. The emission of 5 77(547) mesons was 
observed in the generated sample and its mass was properly reproduced. The 
width of the vr*^ peak is a = 17.8 ± 1.6 MeV/c^ (13.1%). This value is similar 
to the one obtained for the real data [a = 22.6 ± 2.9 MeV/c^, 16.1%), the 
difference being probably due to the larger background at the earth surface. 
The shape of m^-y distribution presents deviations from the gaussian shape as 
it was found with the real data. 

After the inspection of Monte Carlo events coming from the tails of the two 
photons invariant mass distribution it was found that some vr^s were emitted 
from a secondary neutral vertex close to the primary one. Some of them were 
the products of A or K meson decays. Wrong assignment of the decay vertex 
obviously leads to a shift in the reconstructed mass. In the LAr TPC the good 
spatial resolution allows to distinguish majority of the secondary verteces if 
their distance from the main hadronic interaction vertex is larger than about 
1 cm. The problem is much more severe in water or sandwich-like detectors 
which are characterized by much worse spatial resolution. It should be noted 
that the expected percentage of the all secondary neutral pions was found 
to be as large as about 25% at 2 GeV and 43% at 5 GeV in FLUKA based 
generator. 
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Fig. 16. Two photon invariant mass from events initiated by muon neutrino NC 
interactions. Pions from both primary and secondary verteces are included; one pion 
events (a), all (one and multi-pions) events (b). 
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8.1 Neutral current events as a background to CC interactions 

An event with 7r° production in NC i/^ interaction can be misidentified with 
an electron from CC interaction if one of the two following configurations 
occurs: 

1) one photon from the decay can be lost when its energy it is too small to allow 
its distinction in the crowded environment surrounding it (E < 20 MeV); 
the calculation based on the decay kinematics shows that the percentage of 
events with such decay photon, being a function of the parent pion energy, 
never exceeds 2.5%. The showers of energy larger than 20 MeV were always 
properly recognized in the simulated events. 

2) The two photon showers overlap: the angle between the two photon showers 
from the decay is sufficiently large to allow to distinguish the two showers 
only for a 7r° energy up to about 1 GeV. The angle decreases from 70.1° 
at 100 MeV to 13.7° at 1 GeV. If one does not require a full 7r° mass 
reconstruction as a pion selection tool the shower overlap is less critical. 
The possible overlapping effect for more energetic pion decays was already 
studied: it was shown [12] that the ^ method can be efficiently applied, 
in LAr detector, to distinguish the beginning of the electron track from the 
beginning of the lepton pair that comes from the photon conversion. 

Moreover, internal conversion due to the Dalitz decay (1.2% of cases) or 
conversion close to the primary vertex (5.5% of cases within the first 1 cm) 
gives raise to misidentification of a photon shower as an electron induced 
cascade in the 6.7% of cases. In the most frequently populated interval of 7r° 
energy (200-=-400 MeV), in about 2% of events the low energy photon can be 
lost. By multiplying this probability by 0.067 it follows that, in our detector, 
the total probability of the 7r° background is expected to be at the level of 

lo-l 



9 Conclusions 

A sample of 212 events with at least two photon showers pointing at the 
apparent interaction vertex were selected from the data collected during the 
surface test run of the ICARUS detector. Methods of measuring the shower 
energy and shower direction were developed and the invariant mass of the 
photon pairs was reconstructed. The position of the maximum in the invariant 
(7,7) mass distribution for the 196 7r° candidates, obtained after a fiducial 
volume cut, was found to be in agreement with the value of the 7r° mass, thus 
confirming the common assumption that the tt'^s are the main source of the 
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photons. The RMS uncertainty of the reconstructed 7r° invariant mass was 
found to be equal to 36.7 ± 3.4 MeV/c^ (27.3%). 

An improved analysis was carried out in order to clean the initial event sample 
from the events measured in a crowded environment caused by other tracks 
and interactions surrounding the 7r° candidates. As a result, the position of 
the maximum in the invariant (7,7) mass for the 97 surviving tt'^ candidates 
was found again to be in agreement with the 7r° mass. Moreover, the RMS 
uncertainty of the reconstructed 7r° invariant mass was found to be equal 
to (22.6 ± 2.9) MeV/c^ (16.1%), which is significantly narrower than the 
mass distribution for the initial event sample. The 7r° mass distribution was in 
good agreement with the distributions obtained for the Monte Carlo generated 
events with neutral pion production in hadronic and neutrino interactions. 
The agreement between the data and the Monte Carlo confirms the validity 
of the shower energy and shower direction reconstruction methods. It also 
demonstrates the potential of the ICARUS LAr TPC, because simple quahty 
cuts were sufiicient to obtain good results in very difficult conditions of the 
surface test, detector commissioning and tuning of the electronics. 
It was found that in the photon energy range from 50 to 5000 MeV the 
electromagnetic shower energy resolution a{E)/E can be approximated by 
the for mula a {E)/E = PXjsfE + P2 with the parameters PI = (0.33 ± 
0.03) ^MeV and P2 = 0.012 ± 0.002. 

To simulate the future measurements in the CNGS beam, neutrino interactions 
at incident energies between 1 and 5 GeV were generated using the FLUKA 
package. Scanning and measurements of the showers pointing to the primary 
and secondary verteces were carried out, and the two photon mass was 
reconstructed for events with different number of showers. It was demonstrated 
that in most cases a proper reconstruction of up to the three pions coming 
from the same vertex was possible. A few events with ?7(547) meson decay into 
two photons were found in simulated events and properly reconstructed. 
The asymmetric decays of 7r°s as a background to Vf. charge current events 
were discussed. It was found that it is likely to loose the photon of energy 
smaller than about 20 MeV but the percentage of events with such a photon 
never exceeds 2.5%. The frequency of events that could be misidentified with 
events containing single electron shower, because of the low energy photon lost 
for either the Dalitz decay or the events with the second photon converting 
close to the neutrino interaction vertex was estimated to be at the level of one 
per mil. 

In summary: our analysis implies very good prospects for the recognition and 
energy measurement of electromagnetic showers during the ICARUS running 
in the Gran Sasso laboratory, where the background due to other signals in 
the neighborhood of neutrino interactions will be negligible. 
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